Simulated ion-sputtering and Auger electron spectroscopy depth profiling study of 

intermixing in Cu/Co 
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The ion-bombardment induced evolution of intermixing is studied by molecular dynamics simu- 
lations and by Auger electron spectroscopy depth profiling analysis (AESD) in Cu/Co multilayer. 
It has been shown that from AESD we can derive the low-energy mixing rate and which can be 
compared with the simulated values obtained by molecular dynamics (MD) simulations. The overall 
agreement is reasonably good hence MD can hopefully be used to estimate the rate of intermixing 
in various interface systems. 



Low-energy ion beams are commonly used in surface 
analysis and for film growth [1,2]. The use of ion- 
sputtering in the controllable production of nanostruc- 
tures and self-assembled nanoppaterns have also become 
one of the most important fields in materials science 
[2-7]. 

The future of nanotechnology ultimately rests on the 
controllable fabrication, integration, and stability of 
nanoscale devices. However, the understanding of the 
fundamental phenomena leading to the formation, sta- 
bility, and morphological evolution of nanoscale features 
is lacking [8]. As the dimensions of the surface features 
is reduced to the nanoscale, many classical macroscopic 
(continuum and mesoscopic) models for morphological 
evolution lose their validity. Therefore the understanding 
of the driving forces and laws governing mass transport 
involved in the synthesis and organisation of nanoscale 
features in solid state materials is inevitable. Unfortu- 
nately the fundamental understanding and the nanoscale 
control of interdiffusion is not available yet [9,10]. In or- 
der to get more insights in the atomic relocation processes 
during postgrowth low-energy ion-sputtering, it is impor- 
tant to measure and to calculate accurately the rate of 
intermixing at the interface. 

In this Letter we will show that using a newly devel- 
oped code for simulated ion-sputtering based on molecu- 
lar dynamics we are able to get mixing rates. We will also 
show that it is possible to extract the mixing rate data 
from AESD. It turns out that the agreement of the two 
methods is reasonably good. In this way we also could 
check the reliability of molecular dynamics simulations. 

The Auger electron spectroscopy depth profiling anal- 
ysis (AESD) [11,12] measurements were carried out on 
a Cu/Co multilayer system. The sample was made by 
sputter deposition on polished single-crystal (111) sili- 
con substrates in a plasma beam sputter deposition sys- 
tem. It was characterized by XTEM and RBS and flat 
interfaces have been found [12]. AESD depth profiling 
was carried out using a dedicated device [13] by apply- 
ing Ar + ions of energy of 1 keV and angle of incidence 
of 10° with respect to the surface of the crystal (grazing 
angle of incidence) . The sample was rotated during ion 
bombardment . 



The depth profiles were measured as a function of the 
sputtering time keeping the bombarding ion current con- 
stant. A STAIB DESA 100 pre-retarded CM A with fixed 
energy resolution was used to record the AES spectra. 
The following AES peaks were detected Cu (60 eV), Cu 
(920 eV) and Co (656 eV, to avoid overlapping). A part 
of a typical depth profile is shown in Fig 1. For clarity 
only the copper is shown in Fig 1. To calculate the con- 
centration the intensity of the copper Auger current was 
normalized to that of the pure copper. The depth scale 
was calculated from the known thickness of the sample 
[12]- _ 

It is generally supposed that during ion mixing the 
atomic movements are similar to that of the usual diffu- 
sion and thus the occurring broadening can be described 
by the same equations [14] . Accordingly in case of a bi- 
layer system the concentration distribution formed due 
to ion bombardment of the interface by a given fluencc 
can be described by the erf function. The variance, a 2 , 
of the erf function determines the extent of ion mixing. 
In many cases a 2 is linearly dependent on the bombard- 
ing ion fluence <E>. In these cases the cr 2 /<!> ratio called 
as mixing rate characterizes the mixing for a given ion 
bombarding condition. Important advantage of using the 
term mixing rate is that it can be in principle directly de- 
rived from the experiments. 

In the case of AESD we measure the Auger current 
of the elements present. The measured Auger current 
depends on the in-depth distribution as follows: Ij — 
+ ij(2)nj + ij('3)K 2 + where Ij is the measurable 
Auger current of element j, ij(k) is the Auger current 
of element j emitted by the atomic layer k, and Kj gives 
the attenuation of the Auger current crossing an atomic 
plane. ij(k) » Xj(k)N(k), where Xj(k) is the concen- 
tration of element j in layer k, while N(k) is the number 
of atoms of the fc-th atomic plane. It is evident that in 
general from a single measured Ij one cannot determine 
the ij(k) values. On the other hand this equation can 
be used to simulate the measured Auger current during 
the depth profiling procedure if we assume an in-depth 
distribution. 

We do not know any experimental measurement of the 
in-depth distributions formed during AESD applying low 



1 




depth 

FIG. 1. The measured depth profile: projectile Ar + , ion 
energy 1 keV, angle of incidence 84° . The concentration of Cu 
is given as a function of the removed layer thickness (measured 
from the interface in A). The diamonds denote the points at 
which a 2 values are determined by AESD. Inset figure: The 
broadening (a 2 ) of the interface as a function of the fluence 
(ion/ A ). The dotted lines show the two extreme of the 
possible slopes (mixing rates). 
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FIG. 2. The simulated filtered (R 2 ) (SAD) as a function 
of the number of the ion impacts (number fluence) at 1 
keV ion energy at grazing angle of incidence using simulated 
ion-sputtering (molecular dynamics). Dotted line is a linear 
fit to the (R 2 ) curves. (R 2 ) is simulated with ion impact 
points at the free surface (9 ML above the interface) and also 
at 4 ML above the interface. 



(0.2 — 2 keV) ion energy. It seems, however, that dy- 
namic TRIM simulation can reliably be used to describe 
AESD [15]. This calculation also provides the in-depth 
distribution during the procedure. It turns out that at 
the beginning of the depth profiling procedure, when the 
interface is still far from the surface the in-depth distri- 
bution can be approximated by the erf function. For the 
evaluation of the experimentally determined depth pro- 
file we will suppose that (i) the interface has an intrinsic 
surface roughness, which can also be approximated by 
an erf function of a . (ii) ion mixing is the only process 
contributing to the broadening of the interface at least 
in the beginning part of the depth profiling, which will 
be studied. Thus for any measured copper Auger current 
we should find xq (the distance of the interface from the 
surface) and a m (the measured variance) values of the erf 
function. Then the variance due to the ion bombardment 
induced mixing is a 2 = (cr 2 ^ — <7 2 ). We derived a 2 at three 
depths (indicated in Fig 1.) to be 16 ±4 A 2 , 72 ± 16 A 2 , 

and 135 ± 35 A . To proceed we must know the num- 
ber of ions causing the broadening at the interface. In 
our experimental arrangement we cannot measure the ion 
fluence. On the other hand we can measure accurately 
the removed layer thickness. Taking the sputtering yield 
from the literature to be Y « 1.2 [12], we can derive the 
curve a 2 vs. fluence $ which is shown in inset Fig. 1. 
Inset Fig. 1 also shows the two limiting slopes; thus we 
can derive the mixing rate being 250 ± 150 A . 

Classical constant volume molecular dynamics simula- 
tions were also used to simulate the ion-solid interaction 
(ion-sputtering) using the PARCAS code [16]. The com- 
puter animations can be seen in our web page [17]. Fur- 
ther details are given in ref. [16,18,19]. We irradiate the 



bilayer Cu/Co (9 monolayers, (ML) film/substrate) with 
1 keV Ar + ions repeatedly with a time interval of 5-20 ps 
between each of the ion-impacts at 300 K which we find 
sufficiently long time for the termination of interdiffu- 
sion, such as sputtering induced intermixing (ion-beam 
mixing) [19]. The initial velocity direction of the im- 
pinging atom was 10° with respect to the surface of the 
crystal (grazing angle of incidence) to avoid channeling 
directions and to simulate the conditions applied during 
ion-sputtering. 

To describe homo- and heteronuclear interaction of Cu 
and Co, the Levanov's [20] tight-binding potentials are 
used [21]. The cutoff radius r c is taken as the second 
neighbor distance. 

Wc randomly varied the impact position and the az- 
imuth angle (f>. In order to approach the real sputtering 
limit a large number of ion irradiation are employed us- 
ing automatized simulations conducted subsequently to- 
gether with analyzing the history files (movie files) in 
each irradiation steps. In this article we present results 
up to 100 ion irradiation which we find suitable for com- 
paring with low to medium fluence experiments. 100 ions 
are randomly distributed over a 40 x 40 A 2 area. The size 
of the simulation cell is 100 x 100 x 75 A 3 including 62000 
atoms. Chanelling recoils are left to move outside the cell 
and in the next step these energetic and the sputtered 
particles are deleted. 

Unfortunately the realistic simulation of the layer-by- 
layer sputter-removal is beyond the performance of the 
computers available. Just to give a hint of the difficul- 
ties, the simulation of few thousands of ion impacts and 
subsequent relaxations should be treated. To reduce the 
computational demand, we carried out two types of sim- 
ulations. In the first case ions are initialized from the 
surface (9 ML far from the interface). 100 ions corrc- 
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FIG. 3. The crossectional view of the ion-sputtered Cu/Co 
bilayer at 1 keV Ar + ion energy (ions are initialized at -4ML) 
obtained by simulated ion-sputtering. The red circles are the 
incorporated ions. 



sponds to $ w 0.063 ion/A (the removal of - 0.35 ML). 
If the ion is initialized e.g. 5 ML below the free surface 
(-4 ML far from the interface) , than we can simulate high 
dose experiments, when nearly 4 ML is sputter-removed. 
Although, the removal of the upper 4 ML has not been 
done, we expect that this kind of an artificial setup of 
ion-sputtering can simulate the real experiment. This is 
because, we expect that the main source of interfacial 
mixing is the progressive enhancement of energy deposi- 
tion at the interface with respect to the depth position 
of the ion-impacts. 

In Fig 2 the evolution of the sum of the square atomic 
displacements (SAD) along the depth direction of all in- 
termixing atoms (R 2 ) = Y^i ~ z i(t = 0)] 2 , where 
Zi(t) is the depth position of atom i at time t, can be 
followed as a function of the number of ions. Using the 
relation a 2 — (R 2 ) it is possible to calculate the mixing 
rate (k — (i? 2 )/$) from the slope of the fitted straigh 
line on the (R 2 ) vs. $ curve in Fig 2. The contribu- 
tion of intermixed atoms to (R 2 ) is excluded in a given 
layer if the in-layer concentration of them is less then ~ 5 
%. AESD can not measure intermixed atoms which has 
very low concentration and which are below the treshold 
sensitivity of AESD. In the 2nd and 3rd Co layers be- 
low the interface we find less then 5 and 1 % Cu which 
are below the sensitivity of AESD measurements. Ion- 
sputtering in the 1st (-9 ML) and in the 5th (-4 ML) lay- 
ers k-QML ~ 317 A and fc_4jwL ~ 397 A , are obtained, 
respectively. The corresponding fluences are $ w 0.063 

ion/A and $ « 1.26 ion/A (the artificial removal of 4 
ML corresponds to $ « 1.2 ion/ A using the F rj 1.2 
ion/atom), respectively. It must also be noted that Cai 

et al. obtained k ~ 400 A 4 using 1 MeV Si + ions and x- 
rax scattering techniques in Co/Cu multilayer [22]. Since 
the ion mixing occur at the low energy end of the cas- 
cade process, e.g. the high energy collisional cascades 
split into low energy subcascades (see e.g. refs. in [18]), 



an agreement is expected for. We also give the measured 
and calculated mixing efficiencies of £ = k/Fp, where 
Fd is the deposited ion-energy/depth. We get £ ps 10 ±2 

A /eV both by experiment and by simulations. On the 
basis of this value we can characterize low-energy ion- 
sputtering induced intermixing in Cu/Co as a ballistic 
interdiffusion process. 

The crossectional view of the ion-sputtered system can 
be seen in Fig 3 as obtained by MD simulations. The 
interface is only weakly intermixed. The incorporated 
ions are also shown in Fig 3. The simulation can also be 
seen as an animation [17]. 

In this Letter we have presented that the combina- 
tion of atomistic simulations with Auger electron spec- 
troscopy depth profiling might be a new efficient method 
to depth profiling analysis of multilayered materials. 
Also, the reasonably good agreement between experiment 
and simulations provides us the possibility of predicting 
interdiffusion properties for various multilayers for which 
no experimental results are available. 
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